How many species have mass Ml* 
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Within large taxonomic assemblages, the number of species with adult body mass M is charac- 
terized by a broad but asymmetric distribution, with the largest mass being orders of magnitude 
larger than the typical mass. This canonical shape can be explained by cladogenetic diffusion 
that is bounded below by a hard limit on viable species mass and above by extinction risks that 
increase weakly with mass. Here we introduce and analytically solve a simplified cladogenetic dif- 
fusion model. When appropriately parameterized, the diffusion-reaction equation predicts mass 
distributions that are in good agreement with data on 4002 terrestrial mammal from the late 
Quaternary and 8617 extant bird species. Under this model, we show that a specific tradeoff 
between the strength of within-lineage drift toward larger masses (Cope's rule) and the increased 
risk of extinction from increased mass is necessary to produce realistic mass distributions for both 
taxa. We then make several predictions about the evolution of avian species masses. 
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Introduction 

For a large taxonomic group under stable macroevolu- 
tionary conditions, how many species have a mass Ml 
This question has wide implications for the evolution 
and distribution of many other species characteristics 
that correlate strongly with body mass, including life 
span , life history, habitat, metabolism and extinction 
risk (Bennett and Owen"sl|l997t iBrownl . Il995t ICalder lit 
ll984HCardillo et al.l . l2005i ). 

Extant species, including mammals, birds, fish and in- 
sects, exhi bit a canon i cal fo r m of the species mass dis- 
tribution (lAllen et all 120061 : iKozlowski and Gawelczvkl 
120021 : IStanlevl . Il973f) . in which the typical mass is an in- 
termediate value; for example, in mammals the typical 
mass is that of the common Pacific Rat (Rattus exulans, 
40 g). Larger or smaller species, in turn, are significantly 
less common, but asymmetrically so: the largest species, 
such as the extinct Imperial Mammoth (Mammuthus im- 
perator, 10 7 g) for terrestrial mammals, are many orders 
of magnitude larger, while the smallest is only a little 
smaller, e.g., the Remy's Pygmy Shrew (Suncus remyi, 
1.8 g). 

The ubiquity of this distribution of species masses 
(Fig. [Tk,) suggests the existence of a universal evolution- 
ary mechanism. A theoretical explanation of this distri- 
bution may shed light both on the interaction between 
ecolo gical and macroevolutionary processes (IStanlev . 
Il975l ). and on long-term trends in species mass ( Alrov . 



"This manuscript is a pre-print version that has not undergone 
final editing. Please refer to the complete version of record in the 
American Naturalist (2008) at 

http: //www. journals .uchicago . edu/toc/an/current 

T Electronic address: aaronc@santafe.edu 



l2000al lbT). including Cope's rule, the oft-studied notion 
that species mass t ends t o increase with i n a lin eage over 
evolutionary time (|Alrovl . Il998t IStanlevl . I1973TI . Clauset 
and Erwin ( 20081 ) recently showed, by comparing ex- 
tensive computer simulations with empirical data, that 
cladogenetic diffusion in the presence of a taxon-specific 
lower limit on mass M m i„ and extinction risks that grow 
weakly with mass, can explain both the canonical form 
described above, and the precise form of the distribution 
of terrestrial mammal masses. 

Here, we present a simplified three-parameter version 
of the Clauset and Erwin (CE) model for the species mass 
distribution and solve this model in the steady state. 
Comparing the predictions of this simplified model with 
species mass data for the same 4002 terrestr ial mammals 
from the late Quaternary (jSmith et alll2003f ) (henceforth 
denoted "Recent" mammals, which includes species that 
have become extinct during the Holocene), we reproduce 
Clauset and Erwin's results. We then show that the 
model's predictions, when appropriately parameterized, 
are also in g ood agreement with data for 8617 extant 
avian species (|Dunning Jrl |2007[) . 



The Clauset-Erwin Model 

Many models of the variation of species body mass over 
evolu tionary time assume a cladogenetic diffusion pro- 
cess (IKozlowski and Gawelczvkl . l2002t iMcKinnevl . Il990t 



1994t IStanlevl . Il973fl where each descendant 
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iMcSheal . 

species' mass Mn is related to its ancestor's mass Ma 
by a random multiplicative factor A, i.e., Mn = \Ma 
(Fig. [T})) , where A represents the total selective influence 
on the desc endan t species' mass from all sources. Clauset 
and Erwin (2008) studied a family of such diffusion mod- 
els, whose form and parameters could be estimated em- 
pirically. Under their model, the evolution of species 
mass is bounded on the lower end by hard physiologi- 



cal limits, perhap s from metabolic ( Pearsonlll948l 195C ; 
West et aL| . 12002ft or morphological constraints (|Stanlev 



19731 ). and on the upper end by an extin c tion r isk that 
increases weakly with mass ([Liow et al. I, 120081) . Near 
the lower limit, however, within- lineage changes to body 
mass become increasingly biased toward larger masses, 
i.e., (In A) increases as M — > M min . 

Using model parameters estimated from ancestor- 
descendant data on extinct North American terres- 
trial mammals s ince the Cretaceous- Tertiary bound- 
ary (jAlrovl . l2008f ). Clauset and Erwin simulated 60 My 
of mammalian body mass evolution and found that the 
predicted distribution closely matches that of 4002 Re- 
cent terrestrial mammal species. They also found that 
simpler diffusion models, e.g., those that omitted either 
the lower limit, an increased extinction risk from in- 
creased mass, or the increased bias toward larger masses 
as M — > M m in, predicted significantly different distribu- 
tions. Thus, these three mechanisms can explain not only 
the canonical form of the species mass distribution, but 
also its taxon-specific shape. 



A Mathematically Solvable Version 

The CE model, however, remains too complex for 
mathematical analysis, even though it omits many eco- 
logical and microevolutionary processes such as interspe- 
cific competition, predation, and population dynamics. 
By simplifying its assumptions slightly, we can formulate 
a diffusion-reaction model of species body mass evolu- 
tion, which retains most of the key features of the CE 
model and which can be mathematically analyzed. This 
model can then be used to make inferences about body 
mass evolution without resorting to laborious simula- 
tions. 

Let /(cc, t) denote the number of species with logarith- 
mic mass x = lnM at time i, i.e., we transform the CE 
model's multiplicative diffusion process into an additive 
diffusion process on a logarithmic scale. The cladogen- 
esis mechanism leads to two offspring species in which 
x — > x + In A. In the continuum limit, a value (In A) > 
(Cope's rule) corresponds to positive drift velocity in 
the equation of motion for f(x,t). In the CE model 
for terrestrial mammals, the drift velocity was estimated 
from fossil data to increase like M 1 / 4 near the lower 
limit on species mass, a feature that Clauset and Erwin 
found was necessary to accurately predict the number 
of small-bodied mammals. (Whether this quarter-power 
form is related to the quarter-power scali ng commonly 
found elsewhere in the body-size literature ( Savage et all 
2004) remains to be seen.) For mathematical simplicity, 
however, we will ignore this dependence at the expense 
of possibly mis- estimating the number of small-bodied 
species. We also assume that selection pressures on body 
mass are roughly independent, implying that the dis- 
trib ution of changes A i s app roximately lognormal (but 
see (|Clauset and ErwiiJ . l20ol) ). 
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FIG. 1 (a) An illustration of the similarity of the body 
mass distributions for different taxon omic groups, using data 
for 4002 Recent te rrestrial mammals llSmith et af 1 120031 ) and 
8617 extant birds l|Dunning Jr.l . I2007I ). For visual clarity, the 
empiric al distributions hav e been smoothed using a Gaussian 
kernel (jWassermanl . [2006) and rescaled to lay on the same 
abscissa interval, (b) A schematic showing the basic clado- 
genetic diffusion process by which species mass varies over 
evolutionary time. The masses of a descendant Mo is related 
to the mass of its ancestor Ma by a random multiplicative fac- 
tor A, which represents the total selective influence on species 
mass from all sources, i.e., Mo = \Ma- 



The feature that two offspring are produced at each up- 
date step corresponds to a population growth term in the 
equation of motion that is proportional to / itself. The 
extinction probability may also be represented by a loss 
term that is proportional to /. For terrestr ial mammals, 
recent empirical studies (|Liow et all 120081 ) support the 
assumption that the probability per unit time of a species 
becoming extinct p e {x) grows weakly with its mass. The 
CE model uses a simple parameterization of this behav- 
ior: p e {x) — A + Bx where B > 0, which corresponds to 
an extinction probability that grows logarithmically with 
species mass. 

Combining these three elements — diffusion, cladogen- 
esis, and extinction — we may write the continuum equa- 
tion of motion f{x 1 1) for the number of species with mass 



3 



x = In M at time t as 



df v d£ _ D (P£ + 
dt dx dx 2 



A-Bx)f 



(1) 



where v = (In A) is the drift velocity (strength of Cope's 
rule) and the variance D = ( In 2 A ) is the diffusion co- 
efficient. (In physics, Eq. |T]) is called the convection- 
diffusion equation or the Fokker-Planck equation.) This 
equation, however, omits the lower limit on species body 
mass, which we incorporate momentarily. 

The time-dependent equation of motion itself may 
be useful for studying evolutionary trends in species 
body mass, or for making inferences about correlated 
extinction or speciation events. For our purposes, how- 
ever, we are mainly interested in its stationary solution. 
Such a steady-state should exist whenever all species 
within the taxon experience roughly the same set of 
macroevolutionary selective pressures, i.e., under stable 
macroevolutionary conditions. To derive this solution, 
we set the time derivative in Eq. (TTJ) to zero to obtain 
f"-lif'+(ot-fa)f = 0, where /i = v/D, a = (1-A)/D, 
(3 = B/ D, and the prime denotes differentiation with re- 
spect to x. We now eliminate the first derivative term 
by introducing / = e^/ 2 tp to transform the steady-state 
equation to 
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ip = 



This equation can be brought into the form 
of the standard Airy's differential equa- 



tion (Abramowitz and Stcgu: 



ry s a 



V/' - zip = 



(2) 



where we introduce the new variable z — f3 x ^x — 
/3~ 2 / 3 {a — /i 2 /4) and the prime now denotes differenti- 
ation with respect to z. The general solution to Eq. ([2]) 
is ip(z) = c\ Ai(z) + C2 Bi(z), where Ai(z) and Bi(z) are 
the Airy functions. Since there can be no species with 
infinite mass we may set ci — 0. The parameter c\ is 
then determined by the normalization of f(x). 

Thus we can now write the species mass distribution 
as: 



f(x) cx e^/ 2 Ai 



V3, 



(1- 



- 2 / 3 rv - fL 



(3) 



Including now the taxon-specific lower limit x m i„ on 
species mass implies the constraint f{x m \ n ) — for the 
steady-state solution, and allows us to eliminate one pa- 
rameter from Eq. ([3]). Using the fact that the first zero 
of the Airy function is located at zq = —2.3381. . . , which 
we now require to coincide with x = x m { n , gives the con- 
straint 

z ={3 1/3 x min -r 2 / 3 (a-fi 2 /£) , 

which we may solve for a. Inserting this result into 
Eq. © yields 



f(x) cx e^ 2 Ai \^ 3 (x - x min ) + z 



(4) 



as the steady-state solution for the species mass distri- 
bution. 

If the lower limit a: m i n is known, this simple-minded 
model has only two parameters: /i, associated with the 
biased diffusion process and (3, associated with the ex- 
tinction process. A positive bias in the diffusion fi > 
(Cope's rule) has several systematic effects on the dis- 
tribution of species masses: it (i) pushes the left tail of 
the distribution away from the lower limit at x m i n , (ii) 
shifts the modal mass toward slightly larger values, and 
(iii) extends the right tail of the distribution. In con- 
trast, increasing (3 implies that species become extinct 
with greater probability for a given mass M, which con- 
tracts the right tail of the distribution. For a given bias 
(i and number of species n, the parameter (3 also set an 
effective upper limit on the expected maximum observed 
mass within the taxon without invoki ng a hard bounda ry, 
e.g., from biomechanical constraints (|McMahorJ , ll973l) . 



Mammalian Body Mass Evolution 

We now test the predictions of this simple math- 
ematical model using empi rical data f or 400 2 Re- 
cent terrestri al mammals dSmith et al. , 2003) and 
8617 birds (|Dunning Jr.l |2007| ). In the former 



case, we take M n 



1.8 g, the size of the small- 



est known mammal, and we estimate v = 0.109 ± 0.021 
and D = 0.508 ± 0.027 (SE) from Alroy's ancestor- 
descen dant data fo r North American terrestrial mam- 
mals ( Alrovl [2008). Incorporating these values into 
Eq. (0| leaves only (3 unspecified. A strong test of 
this model would estimate (3 from fossil data; how- 
ever, while studie s of extinct i on am ong mammals sug- 
gest that (3 > (|Liow et all |2008), current data does 
not appear to be sufficiently detailed to give a pre- 
cise estimate for m ammals. Instead, following Clauset 
and Erwin (|2008l ). we choose [3 by minimizing the 
tail-weighte d Kolmogorov-Sm irnov (wKS) goodness-of- 
fit statistic (jPress et all Il992f) for the predicted and em- 
pirical distributions: 



wKS = max ■ 



\S{x)~P{x)\ 
^P{x)(l-P(x) 



where S(x) is the empirical distribution function and 
P(x) is the predicted cumulative distribution function. 
Thus, small values of wKS correspond to a model that is 
statistically close to the empirical data everywhere. We 
find that two alternative methods of choosing (3, by nu- 
merically matching the modal masses of the model and 
the empirical data or by matching the expected maxi- 
mum mass of the model with the observed maximum in 
the empirical data, produce similar results. 

For the interested reader, details on the preparation 
of t he em pirical data are discu ssed at length by Al- 
roy ( 1998h and Smith et al. (2003) for mammalian fauna, 
and by Dunning Jr. ( 20071 ) for avian fauna. In general, 
body mass estimates were derived using conventional 
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The resulting fit (Fig. [2^) is in good agreement with 
the empirical data, except for a slight overestimate of 
the number of species with mass near 1 kg, an underesti- 
mate of the number near 300 kg, and a slight misestimate 
of the number of very small-bodied species. The devia- 
tions in the right tail are also seen in the CE model and 
may be due to, e.g., phylogenetically correlated specia- 
tion or extinction events in the recent past. The devia- 
tions in the left tail may be due to our omission of the 
mass-dependence in the drift term \i identified by Clauset 
and Erwin; however, incorporating this behavior into our 
diffusion-reaction model is technically non-trivial. 

Thus, the evolution of mammalian species body masses 
can largely be viewed as a simple diffusion process, char- 
acterized by (i) a slight within-lineage drift toward larger 
masses over evolutionary time (Cope's rule), (ii) a hard 
lower boundary on how small body masses can become, 
and (iii) a very soft constraint on large body masses in the 
form of increased extinction risk. Phrased more concep- 
tually, the left tail of the mammalian body mass distribu- 
tion is mainly controlled by the lower limit on mass, while 
the right tail is the result of an evolutionary tradeoff at 
different timescales: over the short-term, within-lineage 
increases in body mass offer selective advantages such 
as better tolerance of resource fluctuations, bett er ther- 
moregulation, better predator avoidance, etc. ( Brownl . 
ll995HCalder Ilj . ll984l) . while they also increase the long- 
term risk of extinction — a tradeoff previously iden- 
tified in the more specific case of carnivorous mam- 



mals ( Van Valkenburghl 119991 : IVan Valkenburgh et all 
120041 . 



FIG. 2 (a) The solution of the diffusion-reaction model of 
species body mass evolution [Equation Q] with (smoothed) 
empirical data for 4002 terrestrial mammals from the late 
Quaternary. Parameters fj, and a; m j n were estimated from em- 
pirical data (see text) on extinct terrestrial mammals since 
the Cretaceous- Tertiary boundary while f3 was chosen to min- 
imize the tail-weighted Kolmogorov-Smirnov (wKS) distance 
between the model and the empirical distribution, (b) The 
{n, /3)-plane of wKS distance of the specified model and the 
empirical data for terrestrial mammals (lighter values cor- 
respond to smaller wKS). The distinctive groove (blue line) 
demonstrates that a systematic relationship between the bias 
parameter \x and the extinction parameter /3 is necessary to 
produce realistic mass distributions (small wKS). The param- 
eter pair used in (a) is marked. Dashes demarcate the line of 
no within-lineage bias (fi — 0). 



techni ques (for examp le, dentition techniques for mam- 
mals ( Damuthl Il99dh ) . For simplicity, differences due 
to sexual dimorphism, geographic variation, etc. were ig- 
nored or averaged out. Although such differences can be 
critical for smaller studies, given the scale of our data, 
in terms of the number of species studied and the wide 
range of body masses, mild misestimates of body masses 
are unlikely to change our conclusions unless they are 
widespread and systematic. 



Avian Body Mass Evolution 

Unlike mammals, data on most other taxonomic 
groups are generally not sufficient to yield accurate es- 
timates of the param eters fi and (3 (but see Novack- 
Gottshall and Lanier (2008)). The evolutionary history 
of mammalian body masses is relatively clear, in part be- 
cause mammalian fossils are relatively plentiful, are often 
sufficientl y well-preserved that body mass estimate s can 
be made (|Damuthl . Il990t IVan Valkenburghl [l990), and 
the distribution of species body masses during an appar- 
ently stable evolutionary period is known. Avian species, 
however, present an interesting case for study using our 
model: the distribution of extant avian body masses 
(Fig. [5^,) is relatively well characterized ( Dunning Jr.l . 
I2007H and evidence of a minimum species body mass 
M m i n is reasonable ( Pearsonl Il950f h However, the avian 
fossil record may be too sparse to yield accurate estima tes 
of n and (3 (|Fountaine et all 120051 : iHone et all I2008H . 

Even without estimates of /j, and 0, however, the diffu- 
sion model can be used to make quantitative statements 
about the general character of avian body mass evolution. 
To demonstrate this, we consider which combinations of 
the parameters /i, (3 produce "realistic" mass distribu- 
tions, i.e., those with a small distributional distance to 
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mass evolution is biased toward lar ger sizes or not (see 
for ex ample (iBokmal . |2003 Im aure rl Il998t iMaurer et all 
1992)) is that a positive bias is not a necessary condition 
for the evolution of realistic mass distributions. Indirect 
tests of the sign of /Lt, b ased either on the mass dis tri- 
bution within subclades ([McSheal . Il994t IWand . l200lh or 
on changes to minimum an d maximum mas ses within a 
subclade over geologic time (jJablonskil . Il997l ) may be ad- 
equate if confounding hypotheses can be eliminated or 
if an appropriate null model is available. On the other 
hand, the most accurate method remains the direc t anal- 
ysis of large amounts of ancestor-descendant data ( Alrovl . 
ll998ll2000btlNovack-Gottshall and Lanierl . 120081 ). prefer- 
ably when derived from as realistic a phylogeny as pos- 
sible. An alternative approach, however, could use the 
systematic relationship between \i and (3 (Fig. (2b) to in- 
fer fi from an empirical estimate of (3. 

Performing the same analysis with the birds data 
(Fig. \3jp) , we find that realistic avian mass distributions 
are also only produced by particular combinations of pa- 
rameters, and also that the groove of minimum wKS val- 
ues follows a systematic relationship with the same basic 
form as that of terrestrial mammals (FigOs), i.e., /i oc (3 
with 9 w 0.5. Like mammals, this groove passes through 
the point /i = 0, implying that this test cannot rule out 
the possibility that Cope's rule (/z > 0) does not hold for 
the evolution of birds. However, there is an important 
qualitative difference between these two grooves: as a 
function of /3, the birds groove grows considerably slower 
than the mammals groove. This indicates that for a given 
diffusion bias fi, a significantly larger extinction parame- 
ter (3 is required to produce a comparably realistic mass 
distribution. In other words, for a particular body mass 
M, the risk of extinction for an avian species should be 
substantially higher than for a terrestrial mammal species 
of the same size. 

Finally, taking arbitrarily the case of no within-lineage 
bias (li = 0) and using the fitted value of /3, we see that 
the predicted mass distribution for birds is generally in 
good agreement with the empirical data (Fig. [3^i) . The 
deviations in the left and right tails may make interesting 
objects for future study. For instance, the deviations in 
the left tail persist across the groove of minimum wKS 
values, even for unrealistically large values of /i. This 
behavior lends some support to the hypothesis that the 
left tail's particular shape is caused by strong physiolog- 
ical constrai nts on body m ass evolution near the avian 
lower limit (|Stanlevl . Il973lh which are not included in 
our model. Additionally, the predicted right tail may 
fit better, with a slightly smaller estimate of (3, were re- 
cently extinct species included, such as the Giant Moa 
(Dinornis robustus, 250 kg). 
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FIG. 3 (a) The steady-state solution of the diffusion-reaction 
model of species body mass evolution [Equation Q], with 
fi = 0, Mmm — 2 g and j3 chosen to minimize the tail-weighted 
Kolmogorov-Smirnov (wKS) distance between the model and 
the empirical distribution, with (smoothed) empirical data 
for 8617 avian species, (b) The (/Lt, /3)-plane of wKS distances 
(lighter values correspond to smaller wKS), showing a sys- 
tematic relationship between p and f3 for producing realistic 
(small wKS) avian mass distributions. The form of this rela- 
tionship is similar to that of mammals (Fig. [2b) , suggesting 
a similar overall macroevolutionary mechanism. 



the empirical distribution. In particular, we compute the 
wKS distance between the model and the empirical data 
over the (fi, /3)-plane and determine the regions that yield 
the best fits. To illustrate this technique in a better un- 
derstood context, we first apply it to the data on Recent 
terrestrial mammals, disregarding for the moment that 
we have an estimate of /i from fossil data. 

The result of this exercise (Fig. [2b) shows that realis- 
tic mammalian mass distributions can be produced by a 
wide, but not arbitrary, variety of biases li and extinction 
risks (3, including no bias at all, i.e., li = 0. This degen- 
eracy, which forms a groove in the (li, /3)-plane following 
roughly li oc f3 e with 9 w 0.5, suggests that part of the 
difficulty in determining for a particular taxon whether 



Discussion and Conclusions 

With respect to birds, our analysis leads us to make 
several concrete predictions about their evolution, (i) As 
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with terrestrial mammals, avian evolution near the lower 
limit M mm is highly constrained, but in ways that that 
differ to some degree from those faced by mammals near 
their lower limit. This implies that avian species evolving 
near this limit should, more often than not, leave larger 
descendants (In A) > and that the strength of this ten- 
dency should increase strongly as M — > M m i n . (ii) The 
body masses of large avian species (M > 20 g) , like those 
for terrestrial mammals, are constrained mainly by ex- 
tinction risks that increase progressively with body mass. 
And, (iii) avian species with a mass M face a significantly 
greater long-term risk of extinction than do mammalian 
species of the same size. This difference in extinction 
risk is large enough that an avian species with a mass of 
roughly 10 5 g should experience an extinction risk com- 
parable to that of a mammalian species roughly 100 times 
larger (10 7 g). 

From a more conceptual perspective, the similarity of 
the results for extant birds and Recent terrestrial mam- 
mals suggests that their evolutionary histories, in terms 
of the processes that govern the variation of species body 
sizes over evolutionary timescales, are fundamentally the 
same. The similarity of these distributions to those of 
other taxonomic groups suggests that this explanation 
may be universal, although further empirical work is 
necessary to substantiate this hypothes is. Indeed, re - 
cent analyses of fossil data for dinosaurs ( Carrand l2006f ) 
seems to support this view. 

On the other hand, although the agreement of the sim- 
ple diffusion-reaction model given in Eq. Q, when ap- 
propriately parameterized, and the observed mass distri- 
butions of Recent terrestrial mammals and extant birds 
is quite good, the model is obviously incomplete in many 
ways. In particular, the number of small-bodied species 
for both mammals and birds is overestimated by the 
model, which is likely because we omitted the increased 
positive bias toward larger sizes identified by Clauset and 
Erwin. Additionally, speciation and extinction events are 
assumed to be independent, and thus the model cannot 
explain fluctuations in the mass distribution caused by 
phylogenctically correlated extinctions, where entire gen- 
era or niches are wiped out by ecological or environmental 
changes. 

We conclude by noting that the model's good agree- 
ment with data suggests that it may be a useful way to 
establish null-expectations in the study of general trends 
in body mass evolut ion (much like diffusio n models in 
population genetics ( Hartl and Clark! 120071 )) in the ab- 
sence of factors such as inter-specific competition, popu- 
lation dynamics, geography, predation, etc. For instance, 
the fully time-dependent formulation in Eq. ([I]) could be 
used to correctly determine the significance of s tatisti- 
cal tre nds in body masses over evolutionary time ( Alrovl . 
l2000af ). Additionally, given the strong correlations be- 
tween body mass and other species characteristics, this 
model of body mass evolution may provide a way to unify 
certain aspects of ecology and evolution. 
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